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Abstract 

It is generally assumed that on sufficiently large scales the Universe is well-described as a 
homogeneous, isotropic FRW cosmology with a dark energy. Does the formation of nonlinear 
cosmic inhomogeneities produce a significant effect on the average large-scale FLRW dynamics? 
As an answer, we suggest that if the length scale at which homogeneity sets in is much smaller 
than the Hubble length scale, the back-reaction due to averaging over inhomogeneities is 
negligible. This result is supported by more than one approach to study of averaging in 
cosmology. Even if no single approach is sufficiently rigorous and compelling, they are all in 
agreement that the effect of averaging in the real Universe is small. On the other hand, it 
is perhaps fair to say that there is no definitive observational evidence yet that there indeed 
is a homogeneity scale which is much smaller than the Hubble scale, or for that matter, if 
today's Universe is indeed homogeneous on large scales. If the Copernican principle can be 
observationally established to hold, or is theoretically assumed to be valid, this provides strong 
evidence for homogeneity on large scales. However, even this by itself does not say what the 
scale of homogeneity is. If that scale is today comparable to the Hubble radius, only a fully non- 
perturbative analysis can establish or rule out the importance of cosmological back-reaction. 
This brief elementary report summarizes some recent theoretical developments on which the 
above inferences are based. 
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1 Introduction 



The Universe that we see around us is lumpy - it has stars, galaxies, clusters of galaxies, superclus- 
ters, sheets, filaments and voids. We do not precisely know from observations what the size of the 
largest structures is; the size beyond which there are no larger structures. On the other hand the 
early Universe is very well-described as a homogeneous, isotropic FRW cosmology [Big-Bang nu- 
cleosynthesis and the relic CMB are evidence of success] and the present Universe is well-described 
as an FRW cosmology with dark energy. How does one reconcile a universe which is observed to 
be inhomogeneous and anisotropic on smaller scales, with a universe that is assumed to be ho- 
mogeneous and isotropic on large scales? Clearly, some way of averaging the matter distribution 
and the related Einstein equations has to be invoked. What is the right way? The true metric of 
the universe is the one produced by the inhomogeneous matter distribution. On large scales, one 
assumes the average of the true metric to be FLRW, constructs the Einstein tensor for it, and uses 
it on the left hand side of the Friedmann equations wherein the matter content on the right hand 
side is a perfect fluid. Because Einstein equations are nonlinear, the Einstein tensor constructed 
from the average metric tensor will in general not be the same as the average of the Einstein tensor 
of the true metric : 



Here < ... > denotes the averaging operation [whose correct definition, for tensors on a curved 
spacetime, is itself a major challenge]; and < g^u >= g^iulpLRW- The correct averaged Einstein 
equations are of course given by the first pair (the equality) in the above set, whereas in cosmology 
we assume the correct equations to be those given by the second pair (by assuming the inequality 
to actually be an equality). This is obviously done because the latter option is infinitely simpler - it 
is straightforward to write the Einstein tensor for the Robertson- Walker metric, but it is impossible 
to find the true metric of the inhomogeneous universe and then average its corresponding Einstein 
tensor. Since the first and the third terms in ([1]) could differ significantly, we might be working 
with the wrong averaged Einstein equation on cosmological scales. This is the averaging problem: 
how to correctly average Einstein equations, and to find out if the neglected terms (the so-called 
back-reaction) can become important in the late stages of an evolving universe, when nonlinear 
structures such as galaxies and clusters form. In particular, can the back-reaction mimic a dark 
energy, and explain the observed cosmic acceleration? 

The problem of averaging of Einstein equations has a long history, and has recently been reviewed 
in an important article by Ellis [Ij. Important contributions to the study of averaging have been 
made in recent times, amongst others, by Buchert [2], Coley [3], Wald [1], Zalaletdinov [5] and 
their collaborators. Specific applications have been developed by Kolb [B], Marra [7], Rasanen 
[S], Sussman [9], Wiltshire pjOj and others. Much of this work, as well as earlier developments, are 
reviewed by Ellis, and we will not enter into details here, except in the context of specific arguments 
developed here. 

If we want to find out the back-reaction on an FLRW universe, it certainly means we are 
taking an FLRW geometry as given on large scales. It is hence necessary to first know what the 
observational evidence for large scale homogeneity and isotropy is, and what is the length scale at 
which homogeneity sets in. 




(1) 
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2 Evidence for large scale homogeneity 



Careful discussions of this issue have recently been given by Clarkson and Maartens [TT] , Maartens 
^12] and Ellis pj. What we have to say below is a summary from these earler works, and is reported 
here because of its significance for the discussion on averaging in the next section. 

It is well known that homogeneity on spatial surfaces can not be established by direct obser- 
vations, because all our observations are on the past light-cone. Hence tests of homogeneity have 
to exploit the following route : isotropy around us along with the Copernican Principle [CP] im- 
plies homogeneity, and hence a FLRW Universe. So one needs to test for isotropy and for CP 
independently. 

The observational evidence for spacetime isotropy around our world-line can be investigated from 
examining the isotropy of the CMB and the galaxy distribution. For a perfectly isotropic CMB, 
all multipoles of the distribution function higher than the monopole, as well their time derivatives, 
vanish. However without CP one cannot deduce the vanishing of the spatial derivatives of the higher 
multipoles, and hence spacetime isotropy about our world-line cannot be deduced without CP. As 
for baryonic matter (along with certain assumptions for the distribution of CDM and dark energy) 
it can be shown that isotropic distribution of the following four matter observables on the light-cone 
implies an isotropic spacetime geometry : angular diameter distances, galaxy number counts, bulk 
velocities and lensing (details and references to original work can be found in [12] )• As pointed out 
by Maartens, it is not known whether almost-isotropy of observations leads to almost-isotropy of 
spacetime geometry. 

Next, one considers what can be inferred about spatial homogeneity, if one assumes CP, and 
considers the following three cases : isotropic matter distribution, isotropic CMB, almost-isotropic 
CMB. If all fundamental observers measure the same isotropic distribution of the four matter 
observables mentioned above, this implies homogeneity, and the Universe is FLRW. It can be proved 
that exact isotropy of the CMB for all observers also implies an FLRW universe. Almost-isotropy 
of the CMB can be shown, via a non-perturbative analysis, to imply an almost-FLRW universe, 
provided some of the time and spatial derivatives of the multipoles are sufficiently small. 

Thus it is clear that the case for an almost-FLRW universe will be strong if observational 
tests support the Copernican Principle. These tests can be carried out by testing the standard 
consistency relations in FLRW geometry. The FLRW curvature parameter which can be inferred 
from geometric measurements is independent of redshift, and a detection of redshift dependence 
of this parameter will indicate departure from homogeneity. A second test is the time drift of 
cosmological redshift, and a third test is a significant difference between the radial and transverse 
BAO scales. None of these tests have yet been carried out, but their eventual execution will play a 
crucial role in confirming large-scale homogeneity. The CP can also be tested by looking for a large 
thermal or kinetic Sunyaev-Zeldovich effect temeperature distortion of the CMB. Also, a large SZ 
effect induced CMB polarization could indicate a violation of CP and hence of homogeneity. 

As of now, there is no evidence against CP, but neither is there clinching evidence for large-scale 
homogeneity. Also, it is not quite clear at what scale homogeneity sets in. If we assume that there 
is homogeneity, and that too at a scale much less than the Hubble scale, say at around 100 MPc, 
then it can be shown (as discussed next) that the cosmological back-reaction is negligible. And the 
ACDM model is then a good description of the present day Universe. On the other hand if there are 
much larger nonlinear structures in the Universe - their formation can then no longer be described 
perturbatively on an FLRW background, and the back-reaction problem will have to be examined 
afresh. 
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3 Averaging in Cosmology and Calculation of Back- Reaction 



Assuming that the scale at which homogeneity sets in is much smaller than the Hubble scale, we 
give three explanations as to why the back-reaction will be small : (i) a simple argument due 
to Peebles [I3]; (ii) our own work [H], [15], [16], [17] which builds on Zalaletdinov's Macroscopic 
Theory of Gravity [MG] [H], [19], [20], [21]; and (iii) the work of Wald and collaborators [22], [1]. 
Similar results obtained by a few other researchers, which support the present inference, are briefly 
reviewed in Paranjape's thesis [23] . 

3.1 An argument due to Peebles ([13] and references therein) 

For nonlinear structures such as galaxies, the Newtonian gravitational potential is of the order 
of the square of the velocity dispersion [about 300 km/sec], i.e. ~ 10~^. Hence the galaxy 
distribution can be described as a perturbation over an FLRW universe. The metric can be written 
as a perturbed FLRW Universe and the Einstein equations can be split into an evolution equation 
for the background scale-factor and the Poisson equation for the perturbed Newtonian potential 
determined by the density contrast (assumed to be provided by non- interacting dark matter). 

In order to find the effect of averaging on the FLRW equations, spatial averages of Einstein 
equations need to be computed to order 0^, in particular for the dominant term which is proportional 
to V0.V0. When this is done, one finds corrections due to back-reaction in both the Friedmann 
equations - corrections in the form of a kinetic energy coming from the mean square velocity 
dispersion, and the averaged gravitational potential energy determined by the density contrast of 
the formed nonlinear structures. Both these correction terms are of the order of a part in a million, 
and hence much smaller than the magnitude of the observed dark energy. 

The discussion by Peebles is patterned in part on the nice work of Siegel and Fry [21] . It seems 
to us that there is room for improvement in this argument : one should not fix the background, but 
allow for the possibility that as perturbations grow, the background around which back-reaction 
should be calculated may itself be changing, because of feedback from the perturbations. One has to 
ascertain that a runaway process leading to breakdown of perturbation theory does not take place. 
This is the study we attempted by applying Zalaletdinov's averaging theory [Macroscopic Gravity] 
to cosmology. Before we summarize our work on applying MG, it will be useful to briefly review 
Buchert's averaging scheme. We do this because the Buchert approach provides simple averaged 
equations, while being less ambitious than MG. Remarkably, the averaged equations that arise from 
MG are very similar to Buchert's averaging equations, enforcing a certain high degree of reliability 
of both approaches, despite their conceptual differences. 

3.2 Buchert's averaging scheme for a dust spacetime 

For a general spacetime containing irrotational dust, the metric can be written as 

ds^ = -dt^ + hij{x, t)dx'dx^ . (2) 

The expansion tensor B* is given by 0* = {l/2)h'^^hkj where the dot refers to a derivative with 
respect to time t. The traceless symmetric shear tensor is defined as a] = 9* — (9/3)(5] where 9 = 9* 
is the expansion scalar. The Einstein equations can be split into a set of scalar equations and a set 
of vector and traceless tensor equations. The scalar equations are the Hamiltonian constraint fl3a]) 
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and the evolution equation for 6 ([Sb 



(3)7^ + - 2(7^ = IQnGp (3a) 

(=^)7^ + e + = 12716-^ (3b) 

where the dot denotes derivative with respect to time t, '•^^T^ is the Ricci scalar of the 3- dimensional 
hypersurface of constant t and o"^ is the rate of shear defined by cr^ = (l/2)(T]cr^. Eqns. ( !3a|) and 
(l3b|) can be combined to give Raychaudhuri's equation 

e + -e^ + 2a^ + 47rGp = . (4) 
3 

The continuity equation p = —Qp which gives the evolution of p, is consistent with Eqns. fl3a|) . f l3b|) . 
We only consider the scalar equations, since the spatial average of a scalar quantity can be defined 
in a gauge covariant manner within a given foliation of space-time. For the space-time described 
by (E]), the spatial average of a scalar \E'(t,x) over a comoving domain V at time t is defined by 

(^) = ^ /" d^xVh'^ (5) 

where h is the determinant of the 3-metric hij and is the volume of the comoving domain given 
by Vx) = Jjj d^xy/h. 

Spatial averaging is, by definition, not generally covariant. Thus the choice of foliation is relevant, 
and should be motivated on physical grounds. In the context of cosmology, averaging over freely- 
falling observers is a natural choice, especially when one intends to compare the results with standard 
FRW cosmology. Following the definition the following commutation relation then holds [2] 

{n-{^) = {m-m{Q) (6) 

which yields for the expansion scalar 

(6)- - (0) = (6^) - (0)^ (7) 

Introducing the dimensionless scale factor a-D = (Vx>/Vx>j„)^^^ normalized by the volume of the 
domain "D at some initial time tin, we can average the scalar Einstein equations f l3ap . (]3bp and the 
continuity equation to obtain 

(0) = 3^ , (8a) 

^y-8v^G(p) + i(7^)= (8b) 

47rG'(p) = Q^, (8c) 

(p)-= -(0)(P)= -3-(p). (8d) 
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Here (7^), the average of the spatial Ricci scalar ^^^TZ, is a domain dependent spatial constant. The 
'backreaction' Qx> is given by 

Qv^l ((e^) - (6)^) - 2{a') (9) 

and is also a spatial constant. The last equation ( 18d|) simply reflects the fact that the mass contained 
in a comoving domain is constant by construction : the local continuity equation p = — 0p can be 
solved to give py/h = po^/ho where the subscript refers to some arbitrary reference time to- The 
mass Ml) contained in a comoving domain V is then Mx> = j^py/hcPx = j^p^y/h^cPx = constant. 
Hence 

which is precisely what is implied by Eqn. (iSdll . 

This averaging procedure can only be applied for spatial scalars, and hence only a subset of 
the Einstein equations can be smoothed out. As a result it may appear that the outcome of such 
an approach is severely restricted, and essentially incomplete due to the impossibility to analyze 
the full set of equations. However one should note that the cosmological parameters of interest are 
scalars, and the averaging of the exact scalar part of Einstein equations provides the requisite needed 
information. A more general strategy would be to consider the smoothing of tensors, which is beyond 
the scalar approach that certainly provides useful information, albeit not the full information. 

The dynamical equations above can be cast in a form which is immediately comparable with 
the standard FRW equations [2]. Namely, 

'dv 4:TtG 



(peff + 3Peff) (11a) 

Pes (lib) 



with Peff and Peg defined as 



- " 16^ " 16^ ' "16^ + 48^- ^^^^ 



A necessary condition for (11 lap to integrate to ( lllbl) takes the form of the following differential 
equation involving Qx> and (7^) 

Qv + 6—Qv + + 2^(7^) = (13) 

and the criterion to be met in order for the effective scale factor a-D to accelerate, is 

Qv > AnG{p) . (14) 

The Buchert scheme has been applied extensively, and in particular can be used to show that 
there indeed are toy cosmological models which when averaged over inhomogeneities can produce 
an apparent acceleration. However, not all Einstein equations are averaged, and one does not have 
an averaged metric here [which we would like to be the FLRW metric]. Macroscopic Gravity can 
achieve that, while reproducing modified Friedmann equtions analogous to the Buchert equations, 
when applied to cosmology. 
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4 Macroscopic Gravity 



This theory is developed comprehensively in the works of Zalaletdinov; briefly introduced in [H], 
and reviewed in Paranjape's thesis [23]. For detailed discussions and the primary interpretation of 
MG, the reader is referred to Zalaletdinov's original papers cited in this article. 

For the purpose of averaging of tensors the key new element which is introduced is a bivector 
{x', x) which transforms as a vector at event x' and as a co-vector at event x. The bivector is 
used to define the "bilocal extension" of a general tensorial object 

P"(x',x) = Wl,{x,x')P'''{x') (15) 

The "average" of P"'{x) over a 4-dimensional spacetime region S with a supporting point x is 

1 



P^{x) = (P") = — / d^x'^-g'P-{x',x) (16) 

\ / ST Jy, 

and this averaging operation preserves tensorial properties. 

There is a certain degree of non-uniqueness in te choice of the coordination bi- vector - the freedom 
coming from the presence of undetermined structure constants in the commutation relations for a 
vector basis in terms of which one can solve for the cordination bivector. The simplest choice is to 
set these structure constants to zero. When that is done, then in a volume preserving coordinate 
system 0"^, [VPC], i.e. one with g{4>"^) = constant, the coordination bivector takes its most simple 
form, namely 

>V/(x',a;) |p,oper=5/. (17) 

The effect of this non-uniqueness on the physical results for averaging in cosmology remains to be 
estimated. Nonetheless, it is noteworthy that the averaged Friedmann equations to be derived from 
this approach are similar to Buchert's and the physical results about the magnitude of the back- 
reaction is identical to the one due to Peebles. This gives confidence in the robustness of the results 
obtained, even though there is freedom in the choice of the coordination bivector. It is also useful 
to note that this bi-vector is different from the Synge bi-tensor which leaves the metric invariant 
upon avaeraging, and hence cannot really be used to average an inhomogeneous geometry. 

Averaged Geometry : the key idea of Macroscopic Gravity is that the average connection 17"^ (x) 

n\ ^ , (18) 

is defined as the connection 1-form on a new, averaged manifold M. Next one defines a correlation 
2-form 

zr, = {n\An'^)^^-n\An^^. (19) 

Denoting R"^ = ( r^f,)^^- , where is the curvature 2-form of the inhomogeneous geometry, and the 
curvature 2-form on the averaged manifold as M"^ can be shown to give 

M", = R\ - T.\\ . (20) 

The inhomogeneous Einstein equations 

- \5tg^^n, = -/.t^^^^) , (21) 
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average out to 



(22) 



CI = (n-, - l^tZl.m) G''- (23) 

G^^ is the metric on the averaged geometry. The correlation 2-form is assumed to satisfy cer- 
tain differential conditions which amount to closure conditions for the above system of averaged 
equations. 

These averaged equations of Macroscopic Gravity carry, in a covariant and non-perturbative 
manner, information about the effect of the underlying inhomogeneities on the averaged geometry. 

5 Application of Macroscopic Gravity to Cosmology 

In order to apply MG to cosmology we start with the assumption that Einstein's equations are to 
be imposed on length scales where stars are pointlike objects (we denote such a scale as I/inhom)- 
The averaging we perform will be directly at a length scale Lflrw larger than about 100/i~^Mpc 
or so. This averaging scale is assumed to satisfy I/inhom ^ -^flrw ^ -^Hubble where i^Hubbie is the 
length scale of the observable universe. The averaging will be assumed to yield a geometry which 
has homogeneous and isotropic spatial sections. In other words, we will assume that the averaged 
manifold 7W admits a preferred, hypersurface-orthogonal unit timelike vector field v"', which defines 
3-dimensional spacelike hypersurfaces of constant curvature, and that v"" is tangent to the trajecto- 
ries of observers who see an isotropic Gosmic Background Radiation. (These "observers" are defined 
in the averaged manifold - we will clarify below what they correspond to in the inhomogeneous 
manifold.) Throughout the rest of this article, for simplicity, we will work with the special case 
where the spatial sections on 7W defined by v"" are fiat. (In principle the entire calculation can be 
repeated for non-fiat spatial sections as well.) One can then choose coordinates (i, x"^), A = 1,2, 3, 
on M. such that the spatial fine element takes the form 

^'^^dsl,,,,, = a'{t)dABdx^dx\ (24) 

where Sab — I ior A — B, and otherwise, and we have v"- — {v*, 0, 0, 0) so that the spatial 
coordinates are comoving with the preferred observers. The vector field v"' also defines a proper 
time (the cosmic time) r such that dr = v"'da = v^df. We will further assume that the averaged 
energy-momentum tensor can be written in the form of a perfect fiuid, as 

r« = p^;%+p<, (25) 

where the projection operator is defined as 

n^.^St + ^v,, (26) 

and p and p are the homogeneous energy density and pressure respectively, as measured by observers 
moving on trajectories (in M) with the tangent vector field v", 

p = nv'va ; P=lnlTl. (27) 
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p and p are observationally relevant quantities, since all measurements of the matter energy density, 
especially those from studies of Large Scale Structure, interpret observations in the context of the 
averaged geometry. An important consequence of the above assumptions is that the correlation 
tensor C^, when expressed in terms of the natural coordinates adapted to the spatial sections 
defined by the vector field v", is spatially homogeneous. This is clear when the modified Einstein 
equations are written in these natural coordinates. 

The existence of the vector field v"" with the attendant assumptions described above, allows us 
to separate out the nontrivial components of the (FLRW) Einstein tensor on in a coordinate 
independent fashion - the Einstein tensor can be written as 

Ux) ^ Etv% ; Ux) ^ ^ {nlEt) , (28) 

where ji{x) and j2{x) are scalar functions whose form depends upon the coordinates used. The 
remaining components given by TT^E^Va and the traceless part of nl^nlE'^, vanish identically. Since 
the energy-momentum tensor in Eqn. f l25l) also has an identical structure, this structure is 
therefore also imposed on the correlation tensor C^. Namely, nlC^Va and the traceless part of 
7r*7r^C^ must vanish. This is a condition on the underlying inhomogeneous geometry, irrespective 
of the coordinates used on either or and is clearly a consequence of demanding that the 
averaged geometry have the symmetries of the FLRW spacetime. 

This leads us to the crucial question of the choice of gauge for the underlying geometry : namely, 
what choice of spatial sections for the inhomogeneous geometry, will lead to the spatial sections of 
the FLRW metric in the comoving coordinates defined in Eqn. (^^7 Since the matter distribution 
at scale Linhom need not be pressure-free (or, indeed, even of the perfect fluid form), there is clearly 
no natural choice of gauge available, although locally, a synchronous reference frame can always be 
chosen. We note that there must be at least one choice of gauge in which the averaged metric has 
spatial sections in the form f l2^ - this is simply a refinement of the Cosmological Principle, and of 
the Weyl postulate, according to which the Universe is homogeneous and isotropic on large scales, 
and individual galaxies are considered as the "observers" travelling on trajectories with tangent w". 
In the averaging approach, it makes more sense to replace "individual galaxies" with the averaging 
domains considered as physically infinitesimal cells - the "points" of the averaged manifold Ai. This 
is physically reasonable since we know after all, that individual galaxies exhibit peculiar motions, 
undergo mergers and so on. This idea is also more in keeping with the notion that the Universe 
is homogeneous and isotropic only on the largest scales, which are much larger than the scale of 
individual galaxies. 

Consider any 3 + 1 spacetime splitting in the form of a lapse function A^(t,x'^), a shift vector 
N^{t, x'^), and a metric for the 3-geometry hAsit, x"^), so that the line element on Ai can be written 

as 

(^)ds^ = -{N^ - NaN"^) dt^ + 2NBdx^dt + hAsdx^dx^ , (29) 

where Na = hAsN^. At first sight, it might seem reasonable to leave the choice of gauge arbitrary. 
However the analysis is then complicated. On the other hand, if we make the assumption that the 
spatial sections on Ai leading to the spatial metric f l2^ on Ai, are spatial sections in a volume 
preserving gauge, then the correlation terms simplify greatly. This is not surprising since the MG 
formalism is nicely adapted to the choice of volume preserving coordinates. The case when the 
gauge is left unspecified is dealt with in our original papers. 
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To begin our calculation, we perform a coordinate transformation and shift to the gauge wherein 
the new lapse function N is given hy N — 1/ y/h where h is the determinant of the new 3- metric Hab- 
In general, one will now be left with a non-zero shift vector N"^; however, the condition N\/h = 1 
ensures that the coordinates we are now using are volume preserving, since the metric determinant 
is given hj g = —N^h = — 1 = constant. We denote these volume preserving coordinates (VPCs) 
by (t, x) = {i,x^) = {t,x,y,z), and will assume that the spatial coordinates are non-compact. For 
simplicity, we make the added assumption that N-^ = in the inhomogenous geometry, so that 
gil— —N'^ — —1/h and qja — 0. The line element for the inhomogenous manifold M. becomes 

HP 

i^)ds' = -T7=-^ + hAB{i, xjdx^dx"" . (30) 

Note that in this gauge, the average takes on a particularly simple form : for a tensor p^j{x), with 
a spacetime averaging domain given by the "cuboid" E defined by 

E = {{i, X, y, z) I -T/2 <t< T/2, -L/2 <x,y,z< L/2} , (31) 

where T and L are averaging time and length scales respectively, the average is given by 



/ dt' / dx'dy'dz' 

Jt-T/2 J-L/2 



p)(t',x',y',z') 



(32) 



where the limits on the spatial integral are understood to hold for all three spatial coordinates. We 
define the "spatial averaging limit" as the limit T — > (or T ^ ^Hubble) which is interpreted as 
providing a definition of the average on a spatial domain corresponding to a "thin" time slice, the 
averaging operation now being given by 



I r+L/2 



L3 



dx'dy'dz' 

L/2 



piit,x',y',z') 



+ O (TL5i,,,J . (33) 



(Note the time dependence of the integrand.) Henceforth, averaging will refer to spatial averaging, 
and will be denoted by (...), in contrast to the spacetime averaging considered thus far (denoted by 
{...) grj.). The choice of a cube with sides of length L as the spatial averaging domain was arbitrary, 
and is in fact not essential for any of the calculations to follow. In particular, all calculations can 
be performed with a spatial domain of arbitrary shape. We will only use the cube for definiteness 
and simplicity in displaying equations. The significance of introducing a spatial averaging in this 
manner is that the construction of spatial averaging is not isolated from spacetime averaging, but 
is a special limiting case of the latter and is, in fact, still a fully covariant operation. 
For the volume preserving gauge, we have 

GAB^{hAB)=a\t)5AB: (34) 

where a and / are some functions of the time coordinate alone. A few remarks are in order on this 
particular choice of assumptions. Apart from the fact that the spacetime averaging operation takes 
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on its simplest possible form ( 132|) in this gauge and allows a transparent definition of the spatial 
averaging limit, it can also be shown that the assumptions in Eqn. (134|) are sufficient to establish 
the following relations : 

Here the second equality arises from the condition g^^ = G''^ which can be assumed whenever the 
averaged metric is of the FLRW form. The last equality follows on considering the conditions 
(r^^) = (FLRW)pa^ obvious notation, (the basic assumption of the MG averaging scheme). Eqn. 
fl55]) reduces the line element on to the form 

iM)^^2 ^ ^ a\t)5ABdx^dx'' . (36) 

a [t] 

The line element in Eqn. ( 136|) clearly corresponds to the FLRW metric in a volume preserving gauge. 
In other words, the (spatial) average of the inhomogeneous geometry in the volume preserving gauge 
leads to a geometry with homogeneous and isotropic spatial sections, also in a volume preserving 
gauge. Note that the gauge in Eqn. (1361) for the FLRW spacetime differs from the standard 
synchronous and comoving gauge, only by a redefinition of the time coordinate. The vector field 
v"" introduced at the beginning of this section and which defines the FLRW spatial sections, is now 
given by 



V 



(a^ 0,0,0) ; f;, = G,,t;^= (^-^,0,0,0,^ . (37) 

Before proceeding to the calculation of the correlation terms and the averaged Einstein equations, 
we briefly describe why it is important to consider the spatial averaging limit of the MG averaging 
operation. The key idea to emphasize is that an average of the homogeneous and isotropic FLRW 
geometry, should give back the same geometry. Since the FLRW geometry has a preferred set of 
spatial sections, it is important therefore to perform the averaging over these sections. Further, 
since the FLRW metric adapted to its preferred spatial sections depends on the time coordinate, it 
is also essential that the spacetime average should involve a time range that is short compared to 
the scale over which say the scale factor changes significantly. Clearly then, averaging the FLRW 
metric (denoted ^^^^^^^gab) given in Eqn. fl36|) (which is in volume preserving gauge) will strictly 
yield the same metric only in the limit T — )■ 0. Namely, for the cuboid S defined in Eqn. fl3T]) 



= , (38) 

which should be clear from the definition of the metric. The result {^^^^^^'(jab) = ^^^^^^dab in the 
spatial averaging limit can also be shown to hold for the FLRW metric in synchronous gauge, where 
the coordination bivector W| can be easily computed using the transformation from the VPCs 
{i, x"^) to the synchronous coordinates (r, y"^) given by 

The transformation (1391) will also later allow us to write the averaged equations in the synchronous 
gauge for the averaged geometry. 

We now proceed to calculating the correlation 2-form Z^^^^ and thereby the averaged Einstein 
equations. 
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6 The averaged cosmological field equations 



We start by defining (in any gauge with = 0) the expansion tensor B;^ by 



(40) 



where the dot will always refer to a derivative with respect to the VPC time i, and h'^^ is the 
inverse of the 3-metric /iab- (This also gives the symmetric tensor Qab = {i/'^N)hAB, which is the 
negative of the extrinsic curvature tensor.) The traceless symmetric shear tensor and the shear 
scalar are defined by 

ai^ei-ie/3)5i ; = , (41) 

where & = 6^ = (1/A^)9f In y/h is the expansion scalar. 

The connection 1-forms u^j = F^^dx^ can be easily calculated in terms of the expansion tensor, 
for an arbitrary lapse function N. Specializing to the volume preserving gauge (A^ = /i~^/^), the 
bilocal extensions fi'j of the connection 1-forms are trivial and are simply given by 



o^(x^x) 



(42) 



Since Gab = dab, the connection 1-forms fi*^- for the averaged manifold are constructed using the 
FLRW metric in volume preserving gauge given in Eqn. f l36|) . and can also be easily evaluated. 



We can now construct the correlation 2-form Z";,^ and from there the correlation tensor : 



Ci 



^ ijb - T^^fe^ 



tjm 



(43) 



Now, the components of the Einstein tensor E'^ for the averaged spacetime with metric fl36|) are 
given by 



El = 3a^H^ ; E^ = = E 



B 
t ' 



+ 3if2 + H 



(44) 



where the peculiar splitting of terms in the last equation is for later convenience. Recall that the 
overdot denotes a derivative with respect to the VPC time t, not synchronous time. In terms of the 
coordinate independent objects introduced in Eqn. ( 128|) . we have 



J2[X) 



+ 3H^]+ H 



(45) 



From the averaged Einstein equations we next construct the scalar equations which in the stan- 
dard case would correspond to the Friedmann equation and the Raychaudhuri equation. These 
correspond to the Einstein tensor components. 



Elv'va=3i{x) ; TT',El + Elv%a = 332{x)+3i{x 



(46) 
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and are given by 



3a^H^ = {kTI - CI) VaV^ 



(47a) 



6a^ (^^ + = {-kTI + CI) {v,v' + irl) 



-K 



(p + 3p) + 2[Q« + Q(2)+5(2)] 



(47b) 



Here Eqn. f l47ap is the modified Friedmann equation and Eqn. fl47bp the modified Raychaudhuri 
equation (in the volOume preserving gauge on Ai). We have used Eqn. fl27j) . with the overbar on 
p and p reminding us that they are expressed in terms of the nonsynchronous time t, and we have 
defined the correlation terms 



Q 



(1) 



h a" 



cfi 



) = {my 



(48a) 



5« = [((3)ri^(3)rg,) 

-{dA{\^^fh)dB{\nVh)) 



(48b) 



h 



(48c) 



5(2) = a''{^h^''dA{\nVh)dB{\nVh)) 



h 



^S^''{dAilnVh)dBilnVh)). 



(48d) 



In defining Q^^) we have used the relation 9^-9^9^ = {2/3)Q^-2(t^. Q^^) and Q^^^ are correlations 
of the extrinsic curvature, whereas S^^^ and S^"^^ are correlations restricted to the intrinsic 3-geometry 
of the spatial slices of Ai. Since the components of are not explicitly constrained we can treat 
the combinations (l/2)(Q(i) + S^^^) = -Cg and 2(Q(^) + + 5(2)) _ ^qA _ ^o^ independent, 
subject only to the differential constraints which we will come to below. 

As discussed in the beginning of Section |5], the remaining components of C^ must be set to zero, 
giving constraints on the underlying inhomogeneous geometry. In coordinate independent language, 
these constraints read 

nlCiva = = ir'Ctv' ; 



„i „b 



-TT^ [^\Ct) = . 



(49) 
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Eqns. (119|) reduce to the following for our specific choice of volume preserving coordinates, 

C°=0 ; Ct = ■ Ci-hi{Cj) = 0, (50) 

It can be shown that the VPC assumption = h^^^"^ reduces the correlations Q*-^-* and S^'^'^ defined 
in Eqns. f l48cp and ( I48d[) . as well as several terms in the explicit expansion of Eqn. (!50|) . to the 
form 

^(^?oo/^r- ^r-,^) - (/^)(r- ^r-,,) . (5i) 

It can be shown that 



h 



mcT)u)- (52) 



An interesting point is that the VPC assumption N = h further allows us to assume {h"^^ /h) = 
{h'^^){l/h) consistently with the formalism. Using Eqn. (J35l) this gives us 

(53) 

This shows that the correlation terms Q*^^-* and iS*-^-* in fact vanish, 

Q(2) = = 5(2) , (54) 
and leads to some remarkable cancellations in Eqns. ( !50l) . which simplify to give 

[{VhQjB ^'^ri^) - -{VhQjK = , (55a) 

^^^(^ef (^)r:^^) - 5^^(^0f ^^^r^^) = o , (55b) 

= («^5«) . (55c) 

These simplifications are solely a consequence of assuming that the inhomogeneous metric in the 
volume preserving gauge averages out to give the FLRW metric in standard form. In general, these 
simplifications will not occur when the standard FLRW metric arises from an arbitrary choice of 
gauge for the inhomogeneous metric. 

In order to come as close as possible to the standard approach in Cosmology, we will now 
rewrite the scalar equations fH7|) (which are the cosmologically relevant ones) after performing the 
transformation given in Eqn. fl39|) in order to get the FLRW metric to the form 



^M)^s2 ^ _^^2 ^ a\T)6ABdy^dy^ ; a(r) = a(t(r)) . (56) 

Since Eqns. (jU]) are scalar equations, this transformation only has the effect of reexpressing all 
the terms as functions of the synchronous time r. Although the transformation will change the 
explicit form of the coordination bivector , this change involves only the time coordinate, and 
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in the spatial averaging limit there is no difference between averages computed in the VPCs and 
those computed after the time redefinition. This again emphasizes the importance of the spatial 
averaging limit of spacetime averaging, if we are to succeed operationally in explicitly displaying the 
correlations as corrections to the standard cosmological equations. The correlation terms in Eqns. 
( I48l) are therefore still interpreted with respect to the volume preserving gauge, but are treated as 
functions of r. For the scale factor on the other hand, we have 



Further writing 
equations (H7|) become 



a3ff = i^.ff™ ; a-f^ + Si/^Ui^l. (57) 

a dr \a la dr'^ 



p{r) = -pit{T)) ■ p{r) = pitir)) , (58) 



Hknw='-^P-1[Q'''+S^''], (59a) 



1 d^a AnG 



N 



a dr^ 3 3 



{p + Sp) + Iq^'K (59b) 



We emphasize that the quantities Q*-^-* and S^^\ defined in Eqns. f l48ap and fl48bp as correlations 



in the volume preserving gauge, are to be thought of as functions of the synchronous time r, where 
the coordinate r itself was defined after the spatial averaging. Such an identification is justified 
since we are dealing with scalar combinations of these quantities. Note that Q*-^-* and iS*-^^ can be 
treated independently, apart from the constraints imposed by conservation conditions, which we 
turn to next. These conservation conditions can be decomposed into a scalar part and a 3- vector 
part, given respectively by 

a-rv'Ct^, = ■ 4Ct,a = 0. (60) 
In the synchronous gauge f l56|) for the FLRW metric, the scalar equation reads 

(a,Q(i) + 6i/FLRwQ^'^) + {drS^'^ + 2ifpLRw5(^)) = . (61) 

We recall that this equation is a consequence of setting the correlation 3-form and the correlation 
4-form to zero, and it relates the evolution of Q^^^ and S^^\ The 3- vector equation (on imposing 
the first set of conditions in Eqn. fH^ ) simply gives d^C'^ = 0, so that C\ = = constant, which 
also implies that = = constant and hence this equation gives nothing new. (We have used 
the relations Cq = C^, C° = d^C^ and Cq = {l/d^)C^ where denotes the nonsynchronous time 
coordinate t.) 

The cosmological equations fl59]) . along with the constraint equations fl5^ and f lM]) are the key 
results of this section. Subject to the acceptance of the volume preserving gauge on the underlying 
manifold Ai they can in principle be used to study the role of the correction terms resulting from 
spatial averaging. 

6.1 A comparison with the averaging formaUsm of Buchert 

The averaging formalism developed by Buchert is based exclusively on the manifold Ai, and there 
is no analog of the averaged manifold 7W in this scheme. Given an inhomogeneous metric on A4 one 
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takes the trace of the Einstein equations in the inhomogeneous geometry, and carries out a spatial 
averaging of the inhomogeneous scalar equations. 

For ease of comparison, we again recall in brief Buchert's construction, by first writing down 
the averaged equations for the simplest case of pressureless and irrotational inhomogeneous dust. 
The metric can be written in synchronous and comoving gauge as 

ds^ = -dt^ + 6ab(x, t)dx^dx^ . (62) 

The Einstein equations can be split into a set of scalar equations and a set of vector and traceless 
tensor equations. The scalar equations are the Hamiltonian constraint f l63ap and the evolution 
equation for O fl63b|) . 

7^ + - = 167rGp , (63a) 

TZ + dtQ + Q'^ = 12710 p , (63b) 

where TZ is the Ricci scalar of the 3-dimensional hypersurface of constant t, O and cr^ are the 
expansion scalar and the shear scalar defined earlier and p is the inhomogeneous matter density of 
the dust. Note that all quantities in Eqns. ( I63l) generically depend on both position x and time t. 
Eqns. fl63ap and fl63b|) can be combined to give Raychaudhuri's equation 

•9*0 + ^©^ + 2a^ + AuGp = . (64) 



The continuity equation dtp = — 0p which gives the evolution of p, is consistent with Eqns. fl63ap . 
( I63bl) . Only scalar Einstein equations are considered, since the spatial average of a scalar quantity 
can be defined in a gauge covariant manner, within a given foliation of space-time. We return to 
this point below. For the space-time described by fl62|) . the spatial average of a scalar \E'(x, t) over 
a comoving domain V at time t is defined by 

(*)^ = -^ / d'xVb^, (65) 

where b is the determinant of the 3-metric b^s and Vt> is the volume of the comoving domain given 
by Vx> = J^d^xy/b. Spatial averaging is, by definition, not generally covariant. Thus the choice of 
foliation is relevant, and should be motivated on physical grounds. In the context of cosmology, 
averaging over freely-falling observers is a natural choice, especially when one intends to compare the 
results with standard FLRW cosmology. Following the definition f l65|) the following commutation 
relation then holds 

dt i^)^ - {dt^)^ = i^v - {<d)v , (66) 
which yields for the expansion scalar 

dt (0)^ - (9,0)^ = {ex - (0)i • (67) 

Introducing the dimensionless scale factor ax> = (Vp/Vpj)^^^ normalized by the volume of the 
domain V at some initial time tj, we can average the scalar Einstein equations f l63ap . fl63bp and the 
continuity equation to obtain 

dt{p)^= -{Q)^{p)^ ; (0)^ = 3^, (68a) 
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dtayV SttG , , 1 



av J 3 '^'^ 6 



dfav\ 4:nG , , 1 



{p)v-7:iQv + {n)^) , (68b) 



av J 3 ""^^ ■ 3 

Here, the 'kinematical backreaction' Qd is given by 



{p)v+oQv. (68c) 



Qi.-^((0%-(0)p)-2(-% (69) 



and is a spatial constant over the domain V. 

A necessary condition for fl68cp to integrate to (]68bp takes the form of the following differential 
equation involving Qx> and {TZ).^, 

mv + 6^Q^ + dt (7^)^ + 2^ (7^)^ = . (70) 
civ O-V 

The equations above describe the essence of Buchert's averaging formalism, for the dust case. 
We note that the remaining eight Einstein equations for the inhomogeneous geometry, which are 
not scalar equations, are not averaged. These are the five evolution equations for the trace-free part 
of the shear, 

d, [ai) = -Qai -ni + \8j, (a^ - ^9^ + 87rGp) . (71) 

and the three equations relating the spatial variation of the shear and the expansion, 

2 

^B||A=30||B- (72) 

Here, TZ^ is the spatial Ricci tensor and, in Buchert's notation, a || denotes covariant derivative 
with respect to the 3-metric. 

In analogy with the dust case, Buchert's averaging formalism can be applied to the case of a 
perfect fluid, by starting from the metric 

ds"^ = -N^de + hABdx^dx^ . (73) 

The averaged scalar Einstein equations for the scale factor aj) are 

3^ + 4;rG (iV2 (p + 3p)>„ = Qj, + Vv, (74) 

QHl-lQnG{N'p)^ = -Qr^-{N'n)^ ; = ^ , (75) 
where the kinematical backreaction Qjy is given by 

Qv = I mQ)\ - {NQ)l) - 2 (N^a^ , (76) 

and the dynamical backreaction Vv is given by 

Vv = {N'A)^+{edtN)^, (77) 
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where A = V jiu'V iU^) is the 4-divergence of the 4-acceleration of the fluid. Eqn. (175|l follows as 
an Ointegral from Eqn. (1741) if and only if the relation 

- IQnG [dt (N^p)^ + ?>Hv (N^ {p + p)) J = , (78) 

is satisfied. There are also the unaveraged equations (which we do not display here) for the shear, 
analogous to the shear equations flTTj) and fl72|) for dust. 

Buchert's approach is the only other approach, apart from Zalaletdinov's MG, which is capable 
of treating inhomogeneities in a nonperturbative manner, although it is limited to using only scalar 
quantities within a chosen 3 + 1 splitting of spacetime. Buchert takes the trace of the Einstein 
equations in the inhomogeneous geometry, and averages these inhomogeneous scalar equations. In 
the context of Zalaletdinov's MG however, we have used the existence of the vector field v"' in the 
FLRW spacetime to construct scalar equations after averaging the full Einstein equations. As far 
as observations are concerned, it has been noted by Buchert and Carfora that the spatially averaged 
matter density (p)^ defined by Buchert is not the appropriate observationally relevant quantity - 
the "observed" matter density (and pressure) is actually defined in a homogeneous space. Since we 
have done precisely this in Eqn. (127|) . we are directly dealing with the appropriate observationally 
relevant quantity in the MG framework. 

Another important difference between the two approaches is the averaging operation itself. 
Buchert's spatial average, defined for scalar quantities, is given (for some scalar '^{t,x'^)) by fl65|) 
above. On the other hand the averaging operation we have been using (given by Eqn. ( I33l) using 
the volume preserving gauge) is a limit of a spacetime averaging defined using the coordination 
bivector , and is different from the one in Eqn. (1651) . 

Most importantly though, Buchert's averaging scheme by itself does not incorporate the concept 
of an averaged manifold Ai (although the work of Buchert and Carfora [2] does deal with 3-spaces 
of constant curvature). In a recent paper we had argued that Buchert's "effective scale factor" 
ax>it) = (Vx)(i(:)/Vx)(tj„))^/^ must be the scale factor for the metric of the averaged manifold, upto 
some corrections arising due to such effects as calculated by Buchert and Carfora. In the present 
work however, it is clear that such a suggestion is necessarily incomplete due to the presence of Eqns. 
f l55|) constraining the underlying geometry. These constraints are in general nontrivial and hence 
indicate that it is not sufficient to assume that the metric of the inhomogeneous manifold averages 
out to the FLRW form - there are additional conditions which the correlations must satisfy. 

To our understanding, Buchert's averaging formalism is a valid aproach, even though it is based 
on a spatial averaging. A central difference from the MG approach is the issue of closure : not all 
the Einstein equations have been averaged in Buchert's approach, but only the scalar ones. This 
puts a constraint on the allowed solutions considered for the averaged equations: (168|) for the dust 
case, and ( I71|) and (1751) for the fluid case. Solutions to these equations must necessarily be checked 
for consistency with the unaveraged equations for the shear. Further, averaging over successively 
larger scales can bring in additional corrections to the averaged equations, as discussed by Buchert 
and Carfora. Also, if one does not wish to identify Buchert's av with the scale factor in FLRW 
cosmology, one is compelled to develop a whole new set of ideas in order to try and compare theory 
with observation. On the other hand, if one does identify aj) with the scale factor, comparison 
with standard cosmology becomes more convenient, but this brings in additional constraints on 
the underlying inhomogeneous geometry. Thus our conclusion is that the Buchert formalism is a 
correct and tractable averaging scheme, provided all the caveats pointed out in this paragraph are 
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taken care of. Also, when these caveats have been taken care of correctly, the Buchert formalism 
is expected to give the same physical results as the MG approach. We recall that in the covariant 
MG approach also, once a spacetime geometry has been identified for the averaged manifold Ai, a 
gauge must be selected for the geometry on the underlying manifold, in order to explicitly compute 
the correction scalars for comparison with observation. 

The advantage of the MG approach is that it accomplishes in a neat package what the Buchert 
approach, with its attendant caveats, sets out to do. In the MG approach, there are no unaveraged 
shear equations, because the trace of the Einstein equations has been taken after performing the 
averaging on the underlying geometry. Since the averaged geometry is FLRW, the shear is zero by 
definition. There is a natural metric on the averaged manifold by construction, the FLRW metric. 
The correlations satisfy additional constraints, given by Eqns. (!55|) . Thus, once a gauge has been 
chosen and if one can overcome the computational complexity of the averaging operation, the 
cosmological equations derived by us in the MG approach are complete and ready for application, 
without any further caveats. 

In spite of these differences, our equations (155]) and flM]) for the volume preserving gauge are 
strikingly similar to Buchert 's effective FLRW equations and their integrability condition in the dust 
case; and in the case of general A^, the role of Buchert's dynamical backreaction Vv in Eqns. fl7^ and 
f lTHj) is identical to that of our combination of (P^^^ +S^'^^). Concentrating on the volume preserving 
case, the structure of the correlation Q^^^ is identical to Buchert's kinematical backreaction Qxi (or 
Qd in the general case). The correlation S^^'^ appears in place of the averaged 3-Ricci scalar (7?.)^ 
in Buchert's dust equations. This is not unreasonable since Buchert's CR^j, can be thought of as 
(7^)^ = 6/cx)/a|,+ corrections, where Gk-ri/a^ represents the 3-Ricci scalar on the averaged manifold 
which in our case is zero, and hence S^^^ represents the corrections due to averaging. Further, 
these similarities are in spite of the fact that our correlations were defined assuming that a volume 
preserving gauge averages out to the FLRW 3-metric in standard form, whereas Buchert's averaging 
is most naturally adapted to beginning with a synchronous gauge. This remarkable feature, at least 
to our understanding, does not seem to have any deeper meaning - it simply seems to arise from 
the structure of the Einstein equations themselves, together with our assumption Df^Z"^^ = 0. 
In the absence of this latter condition, one would have to consider the correlation 3- and 4-forms 
mentioned earlier, and the structure of the correlation terms and their "conservation" equations 
would be far more complicated. 

An entirely different outlook towards his approach has been emphasized to us by Buchert. 
According to Buchert, the absence of an averaged manifold Ai is not to be thought of as a 'caveat', 
but as a feature deliberately retained 'on purpose'. The actual inhomogeneous Universe is regarded 
by Buchert as the only fundamental entity, and the introduction of an averaged Universe is in 
fact regarded as an unphysical and unnecessary approximation. As we mentioned earlier, this is 
probably the most important difference between MG and Buchert's approach. In the latter, contact 
with observations is to be made by constructing averaged quantities, such as the scalars defined 
earlier in this section, and by introducing the expansion factor ax>- The assertion here is that the 
averaging of geometry, as discussed in MG or in the Renormalization Group approach of Buchert 
and Carfora [2j is not an indispensable step in comparing the inhomogeneous Universe with actual 
observations. The need for averaging of geometry is to be physically separated from simply looking 
at effective properties (such as the constructed scalars) which can be defined for any inhomogeneous 
metric. Averaging of geometry becomes relevant if (i) an observer insists on interpreting the data 
in a FLRW template model, so that (s)he needs a mapping from the actual inhomogeneous slice 
and its average properties to the corresponding properties in this template, or (ii) one desires a 
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mock metric, to sort of have a thermodynamic effective metric to approximate the real one. In this 
context it should perhaps also be mentioned that the importance of a thin time-slice approximation 
of spacetime averaging (as opposed to a strict spatial averaging) has been stressed also by Buchert. 



7 Perturbation theory, structure formation, and backreac- 
tion 

We have in hand the machinery to ask the following question : Is cosmological perturbation theory 
stable against growth of backreaction? The answer must be found iteratively. Assume a background 
with perturbations on it, calculate the back-reaction, feed it in the right hand of the modified 
Priedmann equations to find the new background, and so on : 



a 



(0) 



^ 0(0) ^ c^o) ^ ^(1) ^ ^(1) ^ (79) 



Let the perturbed FLRW metric be 

ds^ = [-{1 + 2(j))dri'^ + 2ujAdx'^dri + ((1 - 2^)7^5 + Xab) dx^dx^] . (80) 

We work with a VPC which has no residual degrees of freedom. Further, this VPC is constructed by 
starting from the conformal Newtonian gauge, and by making a steady coordinate transformation. 
This ensures that all averaged quantities are gauge invariant. We evaluate the correlation scalars 
for a given initial power spectrum - standard CDM. 

For a constant nonevolving potential 0(x), and with a power spectrum 

= Aik/HoT-' , (81) 

the back reaction is 

^~-^(10-'). (82) 

The smallness of backreaction holds also for the exact sCDM model thus demonstrating the stability 
of perturbation theory against the growth of back-reaction. 

This analysis ignores contribution of scales that have become fully nonlinear in matter density 
at late times and it is important to ask if structure formation can significantly modify large sacle 
dynamics. 

We studied backreaction in a toy model of spherical collapse, using the LTB solution. The initial 
density is chosen to be 

{(1 + 5,), r<r, 
(1 - 6^), r^<r <r^ (83) 
1, r > , 

We match the initial velocity and coordinate scaling to the global background solution, by 
requiring 



i?(t,,r) = a,r, (84) 
R{ti,r)^aiHir, (85) 
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For the FLRW background we consider an Einstein-deSitter (EdS) solution with scale factor 
and Hubble parameter given by 

a{t) = {t/uf'^ ■ to = 2/{3Ho), (86) 

H{t) = a/a = 2/(3t) , (87) 

with ^0 denoting the present epoch, fixes the initial time as 

ti^2/(3Ho)af\ (88) 

We use Oj = 10~^, so that the initial conditions are being set around the CMB last scattering epoch. 
The mass function M(r) and curvature function k{r) in this LTB solution are given by 

( 1 + < r < 

GM{r) = -Hyl i + s,((rjrf-l), n < r < (89) 
[ 1 + {Sy/r^) (r^ -rl), r > , 

where we have defined a "critical" radius Tc by the equation 

The significance of Vc is brought out by k{r) : 

= — < 5. ((re/r)' - 1) , r, <r<r, (91) 
I {5./r^) {rl - rl) , r > . 

Since 6^,6^ > 0, we have Vc > r* by definition. The following possibilities arise : 

If rc> Ty, then k{r) > for all r, and every shell will ultimately collapse, including the "void" 
region r* < r < r^. If Tc < r^, then k{r) > for r < Tc and changes sign at r = Tc. Hence, the 
region r* < r < Tc will collapse even though it is underdense, while the region r > will expand 
forever. If r^, = r^, then the "void" exactly compensates for the overdensity, and the universe is 
exactly EdS for r > r„. 

Transforming to the perturbed FLRW form : We want a coordinate transformation (t, r) — (r, f) 
such that the metric in the new coordinates is 

ds^ = -(1 + 2(j))dT^ + a\T){l - 2^) {dr^ + f^dQ^) , (92) 

with at least the conditions 

I |< 1 ; I |< 1, (93) 

being satisfied. Since t is the proper time of each matter shell, the quantity dtf is simply the velocity 
of matter in the (r, f) frame (which is comoving with the Hubble fiow) : 

Br 

v.^, (94) 

is the radial comoving peculiar velocity of the matter shells in the (r, f) frame. We showed that the 
required transformation exists, provided matter peculiar velocities remain small, which is consistent 
with what has een shown by other authors, and is true for the observed Universe. 
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In the cosmological equations derived from Macroscopic Gravity we already have in place the 
formalism for calculating the backreaction when the metric is of the perturbed FLRW form. From 
there it follows that the backreaction is very small, in the nonlinear structure formation regime, 
provided matter peculiar velocities are small. It can be argued that this result is independent of the 
assumption of spherical symmetry in the toy model. The situation could be very different though, 
if there are dominant nonlinear structures in today's Universe, comparable to the Hubble radius. 

7.1 Perturbation theory around a background - the shortwave approx- 
imation 

Green and Wald [22] have recently given an analysis of the growth of metric perturbations, assum- 
ing that the metric is always close to a given background, although matter perturbations can be 
arbitrarily large. No averaging of an underlying spacetime geometry is done, and it is assumed that 
there is a homogeneity length scale at around 100 Mpc, much smaller than the Hubble radius. It is 
shown that if the small-scale motions of matter inhomogeneities are non-relativistic, the deviations 
from the background metric are small, and well-described by Newtonian gravity. This result tallies 
with what has een found by others before, including us. It is further shown that subject to the 
matter satisfying weak energy condition, the effect of small scale inhomogeneities on large scale 
dynamics is to produce an effective trace-free stress energy tensor. One might ask if this traceless 
nature of the correction has to do with no averaging over finite volumes being carried out. 

Thus the assumption of non-relativistic peculiar velocities along with the assumption of a homo- 
geneity scale much smaller than the Hubble radius strongly suggest a negligible effect of small scale 
inhomogeneities on the average large-scale dynamics. The first of these two assumptions is well 
supported by observations. There is no observational evidence against the second assumption, but 
nor is it firmly established by observations. If this assumption is correct, either a small cosmological 
constant, or a modification of general relativity on large scales, is indicated by the observed cosmic 
acceleration. If this assumption turns out to be not correct, the effect of inhomogeneities could be 
significant, and remains an important question for further investigation. 
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The list of references below is far from exhaustive, and references to a large number of the 
original papers on the subject can be found in the review articles cited here. 
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